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We investigate the behaviour of heavy QQ- and QQQ-fiee energies in SU (3) pure gauge theory at 
finite temperature. After stating some general properties of the free energies, we discuss the QQ- 
antitriplet and -sextet channels as well as free energies in different QQQ-channels. We perform 
a detailed finite volume scaling analysis of QQ-fiee energies below and above T c . Moreover, we 
also discuss the relation of the QQ-antitriplet to the QQ-singlet free energies. A comparison of 
the effective running coupling indicates that the r-dependence of the free energies differs by a 
factor of two in agreement with the different Casimir factors. We observe that the baryonic free 
energies can be decomposed into a sum of two particle interactions. 
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1. Introduction 

In this work we discuss the behaviour of heavy (static) diquark and baryonic systems at finite 
temperatures in the quenched approximation. We have calculated static diquark and baryonic free 
energies from Polyakov loop correlation functions. The results of the diquark free energies are 
compared to free energies of quark anti-quark pairs and to the baryonic systems. 

2. QQ Free Energies 

2.1 Observables and Simulation Details 

A system of two quarks can either be in a colour antitriplet (anti-symmetric) or colour sextet 
(symmetric) state. For a static diquark system the free energies in the different colour channels can 
be calculated using correlation functions of the Polyakov loop. For the antitriplet and sextet state 
these are given by ^ 

C&r,T) = ^(TrL(O)TrL(r)) - ~<TrL(0)L(r)) (2.1) 

C 6 qq (r,T) = ^(TrL(O)TrL(r)) + i(TrL(0)L(r)>. (2.2) 

and the corresponding static diquark free energies are obtained by 

F x qq {r,T) = -T\nC x qq (r,T) with *e{3,6}. (2.3) 

The free energies calculated in this way describe the change in free energy of a gluonic system 
containing a diquark pair in a specific colour channel. Both correlation functions are gauge de- 
pendent quantities. Following [Q] we fix to Coulomb gauge. Furthermore they are not Z(3) sym- 
metric and therefore vanish when averaging over the Z(3) sectors. We rotate the configurations 
such that the Polyakov loop (averaged over one configuration) lies in the real Z(3) sector. We 
have performed quenched calculations on lattices with various spatial extent using the tree level 
Symanzik-improved gauge action. We compare our results with the heavy quark-antiquark free 
energies calculated in [J|, ||]. 

2.2 General Properties 

In the following the static diquark free energies are renormalized according to the method 
described in [|3|] by renormalizing the Polyakov loops using the same renormalization constants as 
calculated from the quark-antiquark free energies. 

Figure |I|(left) shows F qq (r, T) / \f& and F qq (r, T) / yfc for four different temperatures below T c . 
We observe clearly F qq (r, T) / \fo to be attractive over the entire distance interval calculated here, 
reaching an almost linear behaviour for large distances. Unlike <2<2-singlet free energies at finite T , 
F qq (r, T)/ \fo retains a T-dependence for small distances. F qq (r, T) / \fo is attractive for distances 
ryfc>\.§, at ryfc « 1.0 it has a minimum and becomes repulsive for even smaller distances, 
remaining T-dependent as well. This qualitative behaviour is already known from QQ-oct&t free 
energies. For large distances F q {r, T) / *fa and F q (r, T) / y/o coincide. 
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Figure 1: Free energies of a static diquark system in the anti triplet (black) and sextet (red) channel below 
(left) and above (right) T c on a 32 3 x 4 lattice. 

In fig. fright) we show the antitriplet, F qq (r,T), and the sextet free energies, F qq (r,T), for five 
temperatures above T c . While F qq (r,T) is attractive over the entire distance range, the <2<2-sextet 
free energies are repulsive for all distances. At large separations the free energies in both colour 
channels tend towards the same asymptotic value. In this limit the colour sources get screened 
independently of their colour orientation. 

2.3 Thermodynamic Limit 



As mentioned in [^J, we rotate the configurations such that the Polyakov loop (averaged over 
one configuration) lies in the real Z(3) sector, which for the Polyakov loop is equivalent to calcu- 
lating (|TrL|). While this expectation value is finite on finite volumes, it vanishes in the thermo- 
dynamic limit. In the strong coupling limit its volume dependence can be understood in terms of 
a random walk [||] that predicts (|TrL|) ~ V~ l l 2 , which we assume to be still valid in the confined 
phase of quenched QCD, and furthermore should hold for C qq (r, T) as well. To check the validity 
of the scaling behaviour of the Polyakov loop and its correlation functions we have performed cal- 
culations on lattices of size N a x N T with fixed temporal extent N t = 4 and various spatial extents 
N a = 16,24,32 and 48. The results for (|TrL|) and C 3 qq (r,T) at a temperature of T /T c = 0.959 
for some values of the separation r are shown in fig. ||(left) and similar for T <T C . They clearly 
show that the expected scaling behaviour holds for all temperatures below T c for the Polyakov loop 
as well as for the correlation functions at all distances. In the thermodynamic limit all non Z(3) 
symmetric observables vanish, despite the fact that we have rotated the gauge field configurations 
to a specific Z(3) sector before calculating this observables. Deviations from a linear behaviour in 
fig. ||(left) are visible for large distances r on the smallest lattice (N a = 16) which may show the 
influence of the periodic boundary conditions. 
We have used the data for a fit of the scaling ansatz 

cl q (r,T)oc Na MrT) oc V - v ^ T \ (2.4) 

The results for the exponent v at temperatures below 0.98 T c shown in fig. fright) are in agreement 
with the expected value of 0.5 within errors. The increase near T c may be explained by an increase 
of the correlation length near the phase transition and indicates that the lattice volumes used may be 
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Figure 2: left: scaling of CL(r,T) for several distances r!T and (|TrL|) with respect to N a foiT/T c =0.959, 



right: fit parameter v(r, T) over T/T c for C qq (r,T) and (|TrL|). 

too small near the critical temperature to observe the correct asymptotic scaling. In the deconfmed 
Phase Z(3) symmetry is spontaneously broken and the expectation values of the Polyakov loop and 
the <2<2-correlation functions are finite in the infinite volume limit. This is shown in fig. |](left) 
where no volume dependence can be identified for all analyzed observables. 



2.4 Comparison to QQ Free Energies 

In perturbation theory the free energy of a quark-antiquark (diquark) pair are in leading order 
given by 

3 



and FUr,T) 



r " ' ' 3 r 

for rT <C 1 and in high temperature perturbation theory {rT 3> 1) 



(2.5) 



4o(rJ) 



e-O and F^rJ) 



2 a(r,T) (T)r 
3 r 



(2.6) 



In order to compare the <2<2-singlet and the <2<2-antitriplet free energies, we compare the ef- 
fective couplings of <2<2-singlet and the <2G- an titriplet, defined by 



a qq (r,T) = --r 



3 2 dF^(r,T) 



dr 



and aUr,T)=--r 



3 *dFl(r,T) 



qq\ 



dr 



(2.7) 



In this way all constant contributions are removed. In fig. fright) we show oc qq (r, T) and oc qq (r, T) 
for temperatures below and above T c . We see, that a^(r,T) and a qq (r,T) agree within errors up 
to large distances. Therefore we conclude, that 2AF qq (r, T) ~ AF qq (r, T), where AF X := F x (r, T) — 
F x (r^°o,T). 



3. QQQ Free Energies 
3.1 Observables 

We have also investigated correlation functions composed of three Polyakov loops represent- 
ing a baryonic system of static quarks on the same configurations as for the diquarks. The colour 
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Figure 3: left: scaling of CL(r, T) for several distances rT and (|TrL|) (scaled by c — 3/8) with respect to 
1 /V for T/T c — 1.682, right: effective couplings T) and a^Jr, T) over r^fa on a 32 3 x 4 lattice. 



state of a QQQ system can be in an anti-symmetric singlet, two symmetrically mixed octet and in a 
totally symmetric decuplet colour state. The observables can be found in [^] and are gauge depen- 
dent quantities, so that we again have fixed to Coulomb gauge. Here we present only calculations 
where the three Polyakov loops are situated at the vertices of an equilateral triangle with the edge 
length R. 



3.2 General Properties 

In fig. gleft) we show the QQQ free energies in the different colour channels and in the average 
channel at 6T C from a 32 3 x 8 lattice. We observe that all channels show screening behaviour 
and approach the same constant for large distances. The singlet free energy is strongly, the octet 
weakly and the average even weaker attractive and the decuplet free energy is strongly repulsive. 
For small distances the singlet free energy coincides with V qqq (T = 0) = 1.5V (T = 0) becoming 
r -independent like the £2<2-singlet free energies. 



3.3 Comparison to QQ Free Energies 

If one assumes that the interaction of the three static quarks of the baryonic system is domi- 
nated by a sum of two particle interactions, the <2<2<2-singlet free energy can be expressed by the 
diquark and quark free energy by 

Fj m (r,T)=3F* t (r,T) + 3F q (T). (3.1) 

In fig. fright) we show the <2<2<2-singlet free energy F qq Jr,T) and three times the <2<2-antitriplet 
free energy F qq (r, T) above T c from a 32 3 x 4 lattice compared to the ansatz (3T ). We observe both 



channels to coincide at all distances and temperatures. We therefore conclude that the free energy 
of the QQQ-^m^&t above T c as a function of distance r can be understood in terms of a sum of two 
particle interactions of QQ pairs being in the colour antitriplet state and the solely T-dependent 
free energy of the static non-interacting quark F q (T). 
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Figure 4: left: Free energies of a static baryonic system in the different colour channels atT/T c = 6 from a 
32 3 x 8 lattice, right: comparison of F qqq (r,T) (red) and 3F qq (r, T) + 3F q (T) (black) above T c from a 32 3 x 4 
lattice. 

4. Conclusion and Outlook 

In this work we discussed the behaviour of heavy (static) diquark and baryonic systems at 
finite temperatures in the quenched approximation. We have calculated static diquark and baryonic 
free energies from Polyakov loop correlation functions in Coulomb gauge. We observed the QQ- 
antitriplet and -sextet free energies to show an behaviour qualitatively similar to that of <2£2-singlet 
and octet free energies. We established C qq (r,T) and (|TrL|) in the thermodynamic limit to vanish 
like 1/yV below T c and stay constant above T c . Furthermore we found the effective couplings 
GC qq -(r, T) and CC qq (r, T) to agree within errors up to large distances. We have shown the free energies 
of the different colour channels of the QQQ system above T c and compared the <2<2G-singlet free 
energy to the <2<2-antitriplet free energy and found, that the free energy of the <2<2<2-singlet above 
T c as a function of distance r can be understood in terms of a sum of two particle interactions of 
QQ pairs being in the colour antitriplet state and the solely T-dependent free energy of the static 
non-interacting quark. 
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